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Abstract. The aim of this paper is to present the analysis of influence of defects in 
ID photonic crystal (PC) on the density of states and simultaneously spontaneous 
emission, in both spatial and frequency domains. In our investigations we use an 
analytic model of ID PC with defects. Our analysis reveals how presence of a defect 
causes a defect mode to appear. We show that a defect in ID PC has local character, 
being negligible in regions of PC situated far from the defected elementary cell. We 
also analyze the effect of multiple defects, which lead to photonic band gap splitting. 
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1. Introduction 

Photonic crystals (PhCs), since first proposed by Yablonovitch ([1]) 
and John ([2]), have attracted considerable attention. These artificial 
structures are materials patterned with a periodicity in dielectric con- 
stant, which can create a range of frequencies (called a photonic band 
gap) for which light is forbidden to exist within the interior of the 
crystal. These photonic band gaps (PBGs) lend themselves to numer- 
ous applications in linear, nonlinear and quantum optics. It has been 
predicted and confirmed experimentally that photonic crystals allow 
to modify and control spontaneous emission of an exited atom due to 
modification of the density of quantum states ([3, 4]). In particular, it 
is well known, that suppression of spontaneous emission is important 
in keeping a system in exited states and can result in the reduction of 
noise in optoelectronic devices. Localized electromagnetic modes can be 
achieved by introduction of defects within the photonic crystal lattice 
([5, 6]). The nature of the defect decides on the shapes and properties 
of localized photonic states in the gap. A point defect could act like 
a microcavity, a line defect like a waveguide. 

In this paper, we study the effect of defects on the density of states 
and simultaneously spontaneous emission rate in one-dimensional pho- 
tonic crystal. We begin by presenting an analytical model of density 
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Figure 1. Structure of a ID photonic crystal. 

of states which is called the effective resonator model (Section 2). In 
general it is well known, that the microcavity affects density of states 
and spontaneous emission rate. In our approach effective cavity model 
is used to study density of states and spontaneous emission rate in finite 
photonic crystal. Moreover, the influence of local defects on density of 
states in spectral and spatial domain is considered. Presented model 
also takes into consideration direction of propagation and therefore 
does not neglect anisotropic nature of photonic crystals. In Section 3 
we briefly derive the analytical expression for mode spectrum which is 
related to the density of states. This relation allow to analyze proper- 
ties of one-dimensional photonic crystals with defects. The influence of 
defects in ID photonic crystal on the density of states and spontaneous 
emission, in spatial and frequency domain is investigated in Section 4. 
Conclusions can be found in Section 5. 



It is an easy task to calculate reflection coefficient of a lattice built of 
dielectric layers, which extend towards infinity in x and y directions 
(where coordinate system is taken as shown in Fig. 1). Following [7], 
in every layer, for each of two polarizations, a solution of Maxwell 
equations can be written as a superposition of two plane waves: 



2. Effective resonator model 
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where i stands for the number of the layer and 
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These plane waves are bound with plane waves in the next layers by 
continuity conditions, which can be written as a matrix equation ([7, 
8,9]): 

Mi M , (3) 
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where Mj is the continuity condition matrix: 
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where /c z is the wave vector's z component in the ith layer, L z is 
the ith layer's width and nu\ is the ith layer's refractive index (more 
precisely, Eq. (3) is a relation between amplitudes of fields in different 
points of space - left edges of layers) . Amplitudes of plane waves outside 
the lattice are therefore related to each other: 
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If ao is an amplitude of the plane wave incident on the lattice, then 6/v+i 
corresponds to the plane wave reflected in the infinity. This reflection 
does not occur, hence bw+i = 0. Therefore: 



where 



and as a result: 
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The reflection coefficient of the whole lattice is then given by 

r= bo__ M 2i 

with the assumption, that amplitudes ao and bo describe incident and 
reflected plane waves in a "virtual" layer with the same refraction coeffi- 
cient as the external medium and width = 0. The above expression 
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Figure 2. Illustration of effective resonator model. 

is purely analytical, though for lattices with diversified layers it would 
take a very complicated form if written as an explicit function of wave 
vector, polarization and lattice parameters, which are hidden in M21 
and M.22- 

Having derived (9), we can calculate reflection coefficient of any 
number of dielectric layers, what allows us to treat any randomly 
selected layer of a lattice as an effective resonator. The idea of this 
model is presented in Fig. 2. Formally, a lattice can be treated as 
two lattices and the layer joined together. It is possible to replace 
the lattices on both sides of the layer with mirrors having the same 
reflection coefficients r^ffe) and rR £ (k), where e denotes polarization. 



3. Mode spectrum 

Let E ex be a plane wave exciting a resonator limited by two parallel 
planes with reflection coefficients and tr (direction of this vector 
determines polarization of the wave). Because of reflections, there ap- 
pears in the resonator an effective plane wave E e ^ . We introduce an 
operator Q, which establishes a relation between them ([10, 11]): 

E eff = QE ex . (10) 
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We call a plane wave E e [kj with polarization e and wave vector k 
a mode of the resonator, if it is an eigenvector of Q: 

QE e (k) =A e (k)E £ (k) . (11) 

Eigenvalue A e (jej is proportional to mode spectrum p e (j?J , which we 
define as: 

„ A«(£) 

*« = • 

Mode spectrum describes distribution of modes in wave vector domain, 
which is equivalent to frequency and direction of propagation. We 
consider the following relation between mode spectrum and density 
of states (DOS): 

dN(k) =Y,k 2 ^ dn Pe (Pj ^ dk . (13) 

In free space no reflections occur and therefore the operator Qf s is 
an identity, with all eigenvalues equal to 1. Hence, mode spectrum is 
constant in free space, and with taken normalization we obtain from 
(13) valid and well known formula for DOS in free space: 

dN fs (k) = %dk . (14) 

Let E e {k^ be a mode of the considered resonator. If the amplitude 
of the excitation is Eq, to find the amplitude E of the effective plane 
wave one has to sum all the reflected waves: 

oo oo 

E = E + Y / Eo (r L r R e 2ik * L *) J +J2 E ° fayr™' 1 ") 3 , (15) 

3=1 3=1 

where L z is the resonator's width. The second sum arises because con- 
tinuity conditions apply on both mirrors, and it can be easily checked, 
that for a plane wave incident on a plane between two different media 
the change of sign of k z (assuming that z axis is perpendicular to 
the plane) causes the reflection coefficient to become its conjugate. 
It follows, that mode spectrum for the resonator is described with the 
formula: 



1 1 - \r L ,r Re \ 2 



8vr 3 l + \r Le r Re \ 2 - 2Re{r Le r Re e 2ik ^ L -} ' 
where reflection coefficients depend on polarization and wave vector. 
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Figure 3. Mode spectrum for perpendicular incidence in a PC with no defects versus 
frequency / and number of elementary cell N . 

4. Defects in ID photonic crystal 

Formula (16) allows to analyze properties of layers of a lattice, in 
a particular case - a ID photonic crystal. In this kind of structure 
there appear photonic band gaps, regions of frequency for which no 
mode is allowed. In terms of mode spectrum a photonic band gap is 
a region, for which mode spectrum is equal to 0. However, in a structure 
built of finite number of layers it is not strictly possible, because the 
structure is not able to fully eliminate coupling with electromagnetic 
field outside of it. Mode spectrum of a perfect photonic crystal in the 
region of a band gap for perpendicular incidence (k = k z e z ) is shown 
in Fig. 3 versus frequency / and the number of elementary cell N in 
which there is the n\ layer, for which the calculation has been made. 
For simplicity, it has been assumed that the crystal is surrounded by 
a uniform medium with refractive index n\. 

As it can be seen in Fig. 3, for elementary cells near the middle of 
the crystal mode spectrum is practically equal to in the gap with 
maximal values at both edges of the gap. Outside the gap mode spec- 
trum fluctuates around the free space value. It means, that the layers 
are resonators with resonant frequencies at the gap edges. For layers 
near the crystals edge there is only a slight difference between mode 
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spectrum of the layer and a uniform medium (which is constant in the 
latter case), because they are resonators with a very low quality and 
do not impose a significant change of an amplitude of a plane wave. 

Introducing a defect of width to a layer results in an appearance of 
a defect mode. The size of the defect determines at which frequency 
does it happen - see Fig. 4 - obviously, because a resonant frequency 
of a resonator depends on its width. However, the defect mode appears 
only in a certain number of layers near the defected one, because elec- 
tromagnetic field with frequency from a band gap has an envelope which 
vanishes in consecutive layers - it has a limited range of penetration, 
hence the defect has no noticeable influence on mode structure of layer 
which is sufficiently distant. This behaviour is depicted in Fig. 5 (values 
on z-axis in this figure are limited). 

Figure 6 shows maximal values of mode spectrum for each frequency 
and in each layer. It can be seen, that the defect introduces much higher 
values of mode spectrum than the edges of the band gap. Modification 
of mode spectrum depends on which layer is defected - a similar graph 
for a defect in the middle of the photonic crystal is shown in Fig. 7. 
In this case the defect mode frequency is the same, but values of mode 
spectrum associated with the defect are even higher. 
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Figure 5. Mode spectrum for perpendicular incidence in a PC with 1 defect versus 
frequency / and number of elementary cell N. 
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Figure 6. Maximal values of mode spectrum in a PC with 1 defect versus (a) 
frequency / and (b) number of elementary cell N for defect in 15 th elementary 
cell. 
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Figure 7. Maximal values of mode spectrum in a PC with 1 defect versus (a) 
frequency / and (b) number of elementary cell N for defect in 30 th elementary 
cell. 
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Figure 9. Maximal values of mode spectrum in a PC with 2 defects versus (a) 
frequency / and (b) number of elementary cell N. 

If a photonic crystal has two defected layers, then there appear two 
regions in which influence of defect modes is significant. Overlap of 
these regions results in two defect modes, with two different frequencies, 
instead of one. This case is shown in Fig. 8 and maximal values of mode 
spectrum are plotted in Fig. 9. Adding more defect leads to further 
splitting of the band gap by a new defect band - see Figs. 10 and 11. 



5. Conclusions 

In this paper we have analyzed the influence of defects in a ID photonic 
crystal on the density of states and spontaneous emission, in both spa- 
cial and frequency domains. We have used the effective resonator model 
and derived an analytical expression for mode spectrum. With the help 
of this relation the presented numerical results revealing properties of 
ID photonic crystal with defects have been obtained. In particular, it 
has been shown that a single defect has local character and its influ- 
ence on mode structure diminishes with the distance from the defect 
position. 
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Figure 10. Mode spectrum for perpendicular incidence in a PC with 5 defects versus 
frequency / and number of elementary cell N. 
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Figure 11. Maximal values of mode spectrum in a PC with 5 defects versus (a) 
frequency / and (b) number of elementary cell N . 
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